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Abstract: The need for precise characterization of dual-phase xenon detectors has grown as 
the technology has matured into a state of high efficacy for rare event searches. The Michigan 
Xenon detector was constructed to study the microphysics of particle interactions in liquid xenon 
across a large energy range in an effort to probe aspects of radiation detection in liquid xenon. We 
report the design and performance of a small 3D position sensitive dual-phase liquid xenon time 
projection chamber with high light yield = 15.2pe/keV at zero field), long elecfron lifetime 
(T > 200 /is), and excellenf energy resolution (a/5 = 1% for 1,333 keV gamma rays in a drifl field 
of 200 V/cm). Liquid xenon lime projeclion chambers wilh such high energy resolution may find 
applicalions nof only in dark maffer direcf defecfion searches, buf also in neufrinoless double befa 
decay experimenfs and ofher applicalions. 
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1. Introduction 

Dual-phase liquid xenon time projection chambers (LXeTPCs) are versatile radiation detection de¬ 
vices capable of both calorimetry and imaging of particle interactions within a condensed xenon 
target. Particularly fruitful applications have been found in dark matter searches i-i , neutrino¬ 
less double beta decay detectors |0, gamma-ray physics experiments [0, medical imaging [0] and 
recently in neutron detection for Homeland Security applications [p. 

While detection of particles in dark matter searches requires experiments to use ever increasing 
xenon target volumes (i.e. multi-ton scale) and to retreat deep underground where background 
rates are low, extensive studies of xenon microphysics and LXeTPC detector characteristics are 
achievable with smaller apparatus on the Earth’s surface. Small LXeTPCs that offer high light 
yield, 3D position sensitivity, and operational flexibility make an excellent platform to carry out 
LXe microphysics studies which probe the effects of quanta production in LXeTPCs, such as signal 
quenching, recombination and anti-correlation, extraction efficiency, energy resolution, light yield, 
charge yield, electric vs. nuclear recoil (ER/NR) discrimination, and more. It is most useful to 
study these effects as a function of deposited energy, interaction type, and with varied TPC electric 
fields. 

LXeTPCs are able to achieve precise interaction vertex reconstruction and good calorimetry 
due to dual signal readout capabilities. During particle interactions in a LXe target, both ionized 
and excited xenon atoms are created. Some of these react with surrounding xenon atoms to form 
short-lived excited dimers (excimers) [0], while some ionization electrons recombine with ionized 
xenon atoms, thus quenching a fraction of the charge signal. This leaves an interaction vertex with 
three populations of quanta: direct excimers (from excited xenon atoms), recombination excimers 
(from ionized xenon atoms), and free electrons. The direct and recombination excimers decay 
rapidly ^(10ns) to produce 178nm photons, called primary scintillation light (SI). An electric 
field is applied across the LXe target to drift ionization electrons away from the interaction site. 
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Figure 1. Photograph of MiX experimental apparatus. At left, the (blue) rack equipment monitors, controls, 
and reads out the necessary components in the rest of the system. At center (red frame), the cryostat houses 
the Mix detector. At right, the (white) gas panel provides storage, handling, and purification of the xenon 
gas. Details of operation reside in the text. 


past a gate grid to the liquid-gas interface, where they are extracted in a strong electric field and 
produce electroluminescence in the gaseous xenon m- The scintillation light produced in this 
process is proportional to the number of ionized electrons and is referred to as S2 [ITTll . A recent 
review of liquid xenon particle detectors gives a detailed account of the detection process [IT^ . 

Detailed study of liquid xenon in radiation interactions is necessary to better understand the 
capabilities of LXe as a target in ionization and scintillation detectors. Previous LXeTPCs have 
been constructed with great success, but more characterization is needed to fully understand the 
nature of quanta detection in LXe and explore the practical limits in energy and position resolu¬ 
tion. The Mix detector was constructed to precisely probe signal production in LXeTPCs, further 
investigating LXe target properties. 

In this paper, the basic features of the MiX apparatus are described in section ^ followed by 
a discussion of the light yield and energy resolution for gamma rays over a large energy range in 
the Mix system in section 0. Finally, in appendix ^ we present a theorem regarding reconstruction 
singularities in detectors with symmetric sensor layout that tends to limit the performance of small 
LXeTPCs with few signal transducers. 


2. Apparatus 

The importance of maintaining a stable operational environment for the MiX detector was con¬ 
sidered from the beginning as an integral part of the design. It was deemed necessary to provide 
and maintain not only adequate pressure and temperature control within the cryostat, but also to 
maintain stable operation of the xenon purification system, the LXe level, the high voltage feed- 
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throughs, the electronics components and the physical space. A photograph of the system is shown 
in Fig. |T]. 

Cryogenics The MiX system consists of a cylindrical inner vessel (6-inch diameter, 12-inch high) 
that is sun'ounded with about 10 layers of superinsulation and placed inside a stainless steel vacuum 
jacket that is evacuated to 10^^ Torr to reduce heat load on the chilled inner volume. The cooling 
power is supplied by a 24 W Joule-Thomson (JT) cooler, model ARS MRCM-150-1. The JT cooler 
is coupled to the inner volume by an indium-sealed copper condensing head that penetrates into the 
gaseous xenon chamber. It pulls heat from the condensing head at a constant rate while a Lakeshore 
temperature control unit (model 340) maintains a constant temperature of the condensation surface 
using a PTIOO platinum resistor temperature sensor and two resistive heating elements placed in 
the condensing head. 

During filling, gaseous xenon (GXe) stored in high-pressure aluminum bottles passes through 
a pressure regulator, condenses on the copper condensing head and falls into a stainless steel funnel. 
The LXe then feeds down to the inner vessel where it cools the detector region. An emergency LN 2 
coil is soldered on the copper condensing head. The coil is permanently attached to a LN 2 dewar 
which can supply cooling power in the event of a cooling failure, and adds cooling power during 
initial cool down. 

Gas system A gas panel, shown in Fig. [I], was assembled to handle the 14 kg of xenon during 
storage and normal operation. Two aluminum bottles are hung from load cells for storage of the 
xenon in gaseous form. One bottle houses all of the xenon while the other acts as a relief/buffer 
volume in case the pressure elevates too high in the inner vessel. Eighteen Swagelok low and high 
pressure pneumatic valves are activated by an electronic relay board. The relay board is controlled 
through a serial-USB interface with a slow control computer to allow remote computer control. 
The valves default to a normally closed configuration in the event that USB connection is reset or 
the power cycled. 

Liquid xenon is drawn from the detector chamber into the heat exchanger [1T311 where it evap¬ 
orates. An MX-808ST-S diaphragm pump pushes the GXe through a hot zirconium SALS PS3- 
MT3-R-1 getter. The mass flow rate is controlled by a Hastings HFC-D-302 flow controller. The 
GXe enters the cryostat again, is cooled in the heat exchanger (93% efficiency [1T311 ), flows into the 
inner vessel condensation region and is liquified again. 

Feed-throughs The ability to communicate with and reliably control the outer vacuum and in¬ 
ner xenon space via both optical and electrical means is crucial for operating the TPC. Custom 
and standard feed-throughs have been used to transport electrical and optical signals between the 
cryostat interior and the external systems. Two standard MPF ten pin feed-throughs with 1.33” CF 
flanges carry low voltage in and out of the detector; one for the cryogenic head located in the outer 
vacuum and another for the inner vessel temperature sensor, the level meter, and the TPC anode. 
Two standard four-connector SHV to vacuum BNC (i.e. non-standard BNC) feed-throughs are used 
to transport the PMT high voltage and signals between the inner vessel and the outside world. The 
SHV-BNC and inner vessel ten pin feed-throughs are mounted to a CF cross and elbow to allow 
for all PMT, temperature, and level sensor cabling to enter the inner vessel through a single 2.75” 
CF flange. 
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Figure 2. Left Panel; Schematic cross section of the MiX TPC. Right Panel; Photograph of the assembled 
detector. Three stainless steel etched hexagonal meshes dehne the drift (cathode - gate; 12 mm) and the 
extraction (gate - anode; 7 mm) regions. Four Hamamatsu R8520-406 PMTs are used on the top and one 
Hamamatsu R11410-MOD is used on the bottom. A copper shielding ring protects the bottom PMT. 

The detector feed-throughs supply negative high voltage ^(lOkV) to the TPC gate and cath¬ 
ode grids via customized commercial ceramic feed-throughs (interlocking teflon sleeves surround 
nickel electrodes coupled to 30kV Kapton jacketed cables with gold connectors). The PMTs are 
fed positive high voltage ^(1 kV) via SHV to vacuum BNC with custom vacuum BNC cable con¬ 
nectors ^ A custom fiber optic feed-through (bare fibre optic core, vacuum epoxy, 1/4” stainless 
tubing, Ultratorr fitting) was constructed to cleanly transport optical photons from a pulsed LED 
into the inner vessel for in situ single photoelectron (pe) calibration of all PMTs. 

Slow Control A slow control system, developed for use in PandaX HTTH, was implemented to pro¬ 
vide monitoring and control of important system variables. Parameters such as pressures, tempera¬ 
tures, setpoints, LXe liquid level, GXe circulation rate, and valve states are logged asynchronously 
during twenty second time intervals while the system is running. Plots of system parameters are 
served to a webpage which can be remotely monitored. An uninterruptible power supply (UPS) 
maintains power to the system if an outage occurs. If the monitored quantities deviate from a pre¬ 
determined range, the slow control system sends alarm messages via email and text message within 
one minute. During emergencies, a liquid nitrogen system is activated to release liquid nitrogen 
from an LN 2 dewar to supply additional cooling power to the condensing head. 

Detector The TPC structure, shown in Fig. is hung from the top flange of the inner vessel 
by four stainless steel spring straps. It consists of four top PMTs (TPMTs) and one bottom PMT 

*It was found that the vacuum side of the standard SHV to vacuum BNC will breakdown at 1 kV—3 kV in ^(1 atma) 
GXe when a standard grounded connector is used. Thus, a machined stainless steel connector with a PTFE core insulator 
was designed to mate to the unjacketed Kap3 coaxial PMT cables. 
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Figure 3. Photograph of MiX detector with all three TPC electrodes in place. Each 50 /rm-thick mesh has 
a central 76 mm circular region with chemically etched hexagonal grid structure that provides a 95% optical 
transparency. The large unetched region extending beyond the TPC sensitive volume ensures that fringing 
electric fields induced by the grounded inner vessel surface do not reach the innermost fiducial region. 

(BPMT) with high detection efficiency in the VUV region (> 25% at 178 nm) [ITEl] , three electrode 
grids (anode, gate, and cathode), one shielding ring for the bottom PMT, and an enveloping virgin 
PTFE diffuse reflector. The 133 mm-diameter hexagonal grids were chemically etched by Great 
Lakes Engineering from 50 /im-thick stainless steel sheet to optimize electric field uniformity. A 
2 mm-pitch hexagonal mesh was formed in the central 76 mm diameter. The wires were smoothed 
as a result of the etching, giving a round wire of 50 /im diameter and an optical transparency of 
95% (see Fig. During operation, the anode grid and shielding ring are connected to ground 
while the gate and cathode grids are supplied with negative high voltage to provide the appropriate 
drift and extraction fields. 

Liquid xenon fills the detector from the bottom until it partially occupies the space between 
the gate and anode (the gap region) which contains both gaseous and liquid xenon. A capacitive 
level meter is used to monitor the height of the liquid-gas interface. The drift region consists 
of the space between the cathode and gate and is the sensitive volume for radiation detection. 
Ionization electrons drift past the gate, are extracted from the liquid, and collected on the anode 
after producing electroluminescence in the gas gap region. Details of the detector dimensions are 
listed in Table IB 

The top PMT array consists of four Hamamatsu R8520-406 tubes. The 1-inch square, low 
background, VUV sensitive tubes are arranged in a symmetric four pixel configuration. Each 
TPMT is affixed to a FR4 positive high voltage base and placed in the top PTFE reflector with only 
the 20 mm by 20 mm active photocathode area exposed to the TPC inner volume. The distribution 
of the S2 signal in the four TPMTs is used to reconstruct the XY position of an event. The bottom 
PMT is a 3-inch diameter Hamamatsu R11410-MOD tube. The active photocathode diameter of 
64 mm is matched to the inner TPC diameter (63.5 mm) to allow for full bottom photocathode 
coverage of the TPC cylinder. A Cirlex positive HV base is attached to the BPMT and supplied 
high voltage via Kap3 cabling with custom gold plated HV connections to the BPMT base. 
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Table 1. Detector parameters. 


Detector Parameter 

Value 

Diameter 

62.5 mm 

TPMT to Anode 

10 mm 

Anode to Gate 

7 mm 

Gate to Cathode 

12 mm 

Cathode to Shield Ring 

17 mm 

Shield Ring to BPMT 

3 mm 

Grid Pitch 

2 mm 

Grid Thickness 

50/im 

Grid Transparency 

95% 


The sensitive volume is designed to be surrounded in near An by either PMT photocathode or 
virgin PTFE reflector to optimize the light yield. Virgin PTFE has been shown to have above 95% 
reflectivity at 178 nm []T9|]. The reflector is machined from PTFE sheet and designed to be only as 
thick as is structurally necessary. The reduction of PTFE mass contained within MiX is done to 
reduce the low Z material content (in an effort to reduce neutron energy loss when neutrons enter 
the detector). 

The PMT gains were measured from the single photoelectron (SPE) signal induced by an 
in-situ pulsed green FED system that is coupled to the bare fiber optic cable and pulsed at low am¬ 
plitude with 10 ns-wide pulses. The fiber optic cable terminates in the top PTFE reflector, allowing 
light pulses to diffuse through the PTFE and into the detector volume. To determine the PMT gains 
at PMT bias voltages where the SPE signal would be poorly resolved from dark noise, the FED 
signal is raised to several hundred photoelectrons and held constant while the PMT voltages are 
lowered incrementally to measure the power law parameters for each PMT. In this way, a low gain 
setting of 5.0 X 10^ could be selected for each individual TPMT. The BPMT was maintained at a 
gain of 1.0 X 10^. 

Electronics Signals induced in the PMTs travel along 50fl SHV cables into custom decouplers. 
The decoupled signals enter a flash ADC (CAEN V1724) that is mated to an optical link (CAEN 
A3818A) in a linux machine. A DRS4 prototype board is used as a digital threshold trigger at 8 mV 
offset. Once an event is triggered, a 30 /is time interval is captured at a sampling rate of 100 MS/s 
for the five PMT channels and stored in a circular buffer. Once the buffer is full, the data are 
transferred over the optical link and written in binary format (typical event rates 10 Hz—300 Hz), 
using the server internal disk as an initial file buffer. Upon completion, the file is automatically 
transferred to network storage for backup and offline analysis. 

TPC Operation The Mix detector was able to maintain linearity over a large energy range by 
careful selection of PMT gains and selective use of top and bottom signals. Since the TPMT 
array sees only a small fraction of the SI scintillation light (< 15%) and has high relative noise 
(due to the low gain selection driven by S2), the SI signal is integrated only from the bottom 
PMT (termed BSl). Additionally, at a gain of 10^, the BPMT saturates on S2 signals beyond 
35 kpe (corresponding to « 25 kpe in the summed TPMT array). This means that the BPMT S2 
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Figure 4. A typical event waveform from a 356keV gamma ray from The inserts show the SI signal 

for the bottom PMT, and the S2 signal for the 4 individual top PMTs. 


signal is only linear for electromagnetic events with energy below 50keV at a 200 V/cm drift field. 
Therefore, the S2 signal was chosen to be derived solely from the summed TPMT array (called 
TS2). This selection allows for a linear response of the detector to better than 5% in both BS1 and 
TS2 for electron recoil events of 3 keV to 1.5 MeV at a drift field of 200 V/cm. The TS2 signal will 
become safurated above 1 MeV if fhe drifl field is 1 kV/cm or greater due fo fhe large high-energy 
TS2 signals resulfing from efficient quenching of recombination. 

A characteristic high energy signal from the MiX detector is shown in Fig. ^ Signals are 
found in the 30 /fs summed PMT output with a pulse threshold of 0.2pe/ns for the results reported 
here. The BSl and TS2 signals are typically well formed and easily identified by widfh using a 
cusfom fhreshold-based pulse finding scheme. Therefore, pulse widfh-above-fhreshold is fhe only 
criterion applied fo defermine whefher a signal is characferized as an SI or S2. In particular, a 
signal is characterized as SI if fhe widfh is above 20ns buf less fhan 750ns, and as S2 if fhe 
widfh is 750 ns or above. The selection process can be adjusted fo include a widfh dependence on 
area if efficiency for proper SI identification of evenfs wifh energies beyond 1.5 MeV is desired. 
However, deferminafion of fhe widfh-dependence-on-energy selecfion was nof sfudied in depfh for 
Ibis analysis, and fhe simple widfh-based selecfion algorifhm is used wifh greaf success. 

The XY position of each evenf is reconsfrucfed from fhe S2 signal disfribufion on fhe lop PMT 
array based on a cenfer of gravity algorifhm, as shown in Fig. ||. Nofe lhal fhe units of this space 
are arbitrary, as they represent an optically warped XY plane. The position resolution is good in 
the center of the detector, as seen from the presence of the pattern imprinted from the hexagonal 
anode mesh. Near the edges, however, the pattern is warped and reconstruction accuracy is poor 
and ultimately degenerate (see Appendix ^ for details). The Z position is determined from the time 
delay between the largest BSl and largest TS2 signal, since the drift time of electrons in LXe is 
known as a function of drift field [HDIl . 
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Figure 5. Reconstructed XY position based on the center of gravity method for gamma rays from *^^Cs. 
The dots observed are at a spacing of 2 mm and correspond to unit cells of the hexagonal anode grid. The 
existence of the dots allows for a fiducial volume selection with high precision in the XY plane. At large 
radius, optical warping and position reconstruction degeneracy is observed. The donut-like region (bounded 
by the two inner, red circles) represents the fiducial region used in the energy resolution analysis, while the 
white outer ring represents the region used in determining the a and j3 values. See text for more details. 

Due to non-uniform light collection and light generation inside the sensitive volume, the BSl 
and TS2 signals depend on event location. Position-dependent signal correction is achieved by 
selecting photo-absorption events from a high statistics ^^Na run. Prohle plots are generated in 
XY, then smoothed and inverted to produce XY corrections, shown in the bottom panels of Fig. ^ 
Additionally, Z prohle plots are generated and then ht with linear functions to obtain the Z depen¬ 
dencies, shown in the top panels of Fig. ^ The two distinct regions in the panels above and below 
the gate electrode rehect the BSl and TS2 response to the higher electric held in the extraction and 
the lower electric held in the drift regions, respectively. In the higher held extraction region recom¬ 
bination is suppressed, leading to a suppression in the primary light signal (BSl), an increase in 
the number of drift electrons and a coiTesponding increase in proportional light (TS2). In the lower 
held drift region, more recombination takes place, leading to an increase in the primary light signal 
(BSl), a decrease in the number of drift electrons and a corresponding decrease in proportional 
light (TS2). 

Corrections for both BSl and TS2 are applied as a function of XY and Z. The BSl and TS2 
corrections are applied in two steps since many more mono-energetic events would be needed to 
build a full three dimensional model. The factorization of corrections in this way is not expected 
to produce a strong systematic effect since the hducial volume is selected to be small and centrally 
located. The corrections help the energy resolution of the detector at the sub-percent level. Two 
hducial volumes were chosen, as displayed in Fig. 0. The larger-radius circular region was used to 
determine the anti-correlation axes while the inner symmetric annular region was used to determine 
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Figure 6. Position dependence of signal production in the detector after selection of the 511 keV photo¬ 
absorption peak of ^^Na. Top left: BSl dependence on Z (i.e. drift time). The vertical solid lines indicate 
the positions of the gate and cathode electrodes. Top right: TS2 dependence on Z. Bottom left: BSl 
dependence on XY. Bottom right: TS2 dependence on XY. The XY and Z dependencies are factored to 
produce corrections for both BS1 and TS2 independently. 


energy resolution. The choice for two distinct fiducial regions was based on the the superior energy 
resolution in the inner region, and the more precise determination of the anti-correlation axes due 
to the higher statistical sample of events in the larger region. The electron lifetime, extracted from 
an exponential fit of the TS2 signal (shown in the top right panel in Fig. ^ was above 200 /f s over 
the course of the two month data taking period. 


3. Results 

Five electromagnetic radioactive sources were used to produce gamma-ray events inside the MiX 
sensitive volume in the energy range of 122keV to 1.33 MeV. Furthermore, ^^^Cf was used to 
produce neutron-activated xenon states of and ^^^"‘Xe at 163.9keV and 236.1 keV, respec¬ 

tively, as listed in Table ^ Xenon activation was performed in situ with two ^^^Cf sources placed 
just outside of the outer vessel, taking advantage of the high density of xenon in the liquid state. 
The combined activity of the ^^^Cf sources was 78.7 /rCi (April 2015), resulting in a neutron flux 
of 3.4 X 10^ neutrons/s. No data were collected during the week-long xenon activation. 

Achieving high energy resolution requires many aspects of the detector to operate predictably. 
The detector must be extremely stable over periods of data collection which may last from a few 
hours to a few days. Any change in the electric fields, photomultiplier tube (PMT) gains, or LXe 
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Figure 7. Correlation of the Top S2 and Bottom S1 signals for all the sources used in this study taken at a 
drift field of 200 V/cm. This set of data was used to determine the a and j3 values. From top left to bottom 
right the plots represent spectra for ^^Co, '^'""Xe and ^^®'”Xe, '^^Ba, ^^Na, '^^Cs, and ^‘’Co. Note that the 
TS2 and BSl ranges increase with energy for each plot but the ratio of TS2 to BSl ranges are held constant 
for all sources. 


level may modify the conditions for events in the detector and worsen the observed energy resolu¬ 
tion. Additionally, a fiducial volume of the detector must be chosen in which the energy resolution 
is good. Events in the dataset must be classified as eifher inside or oufside of fhis fiducial volume 
based on precise XYZ position reconsfrucfion. Finally, small variafions in fhe lighf collection and 
signal generafion across fhe defecfor musf be correcfed wifh an XYZ posifion-dependenf signal 
correcfion. 

Since fhere is af leasf one phofo-absorpfion peak associafed wifh each radioacfive source, fhe 
response of fhe MiX defecfor fo mono-energefic evenfs can be measured for each source. Figure |7| 
shows a 2D representafion of fhe TS2 versus BSl signals recorded for each source af a drifl field of 
200 V/cm. The infrinsic evenf-by-evenf anfi-correlalion of lighf (BSl) and charge (TS2) signals in 
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Figure 8. Field dependence of scintillation (BSl as open triangles) and ionization (TS2 as open circles) yield 
and their sum (filled squares) in LXe for 356 keV electron recoils from *^^Ba. The values are normalized to 
unity at 2,000 V/cm. 


LXe is clearly seen for each source. There are two manifestations (macroscopic and microscopic) 
of anti-correlations to be considered. The mean macroscopic anti-correlation [IT^ in BSl and 
TS2 is seen when sweeping the drift field for a specific gamma line energy (as shown in Fig. 

As fhe driff field is increased, ionizafion collection grows and evenfually asympfofes while lighf 
signal drops complemenfarily. The behavior is due fo ionization elecfrons being collecfed rafher 
fhan forming additional SI recombination signal. The microscopic anfi-correlafion, due fo fhe 
evenf-by-evenf sfochasfic nafure of elecfron recoil energy parfifioning info excifafion and ionizafion 
(wifh fhe effecfs of recombinafion), coupled wifh fhe assumed consfanf work funcfion of LXe 
{W = 13.7 eV) 1^-^ and high defector lighf yield, can be seen for fhe various dafa sefs in 
Fig.0. The consfanf nafure of W (fhe average energy fo produce a scinfillafion phofon or fo liberafe 
an elecfron) across energy is approximafely verified, buf precision measuremenfs are necessary fo 
confirm fhis behavior. 

Alfhough bofh fhe BSl and TS2 signals can be used individually fo defermine fhe energy of 
an inferacfion, if is advanfageous fo combine fhe BSl and TS2 signals for energy reconsfrucfion 
fo achieve an opfimized energy resolution. The combined energy scale Eces is defermined by 
combining fhe fwo signals fo form 


where a (also referred fo as gi, PDF) is fhe average phofon defecfion efficiency and j3 (or g 2 , 
ETE ■ EEE ■ SE) is fhe producf of fhe elecfron exfracfion efficiency (EEE), fhe elecfron fransmis- 
sion efficiency while passing by fhe gafe grid (ETE), and fhe single elecfron S2 gain (SE). This 
combinafion of signals provides fhe optimal energy resolufion af fhe energies sfudied in fhis work. 
Figure ^ (righf panel) shows an example of a combined energy specfrum for ^^Na in a 200 V/cm 
driff field. If corresponds fo fhe projection of evenfs in fhe TS2 vs. BSl specfrum (Fig. ^leff panel) 
along fhe S2-S1 anfi-correlafion axis whose slope is given by fhe ratio of —j8/a. Using fhe large 
XY fiducial region and a 9/^5 > Z > 4/is driff time cuf, an energy resolufion (a/E) of 1.62% for 
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Figure 9. Energy resolution for large fiducial volume selection as shown in Fig. 5, but no XYZ position 
corrections applied. Left panel; Correlation of the Top S2 and Bottom SI signals for ^^Na in a 200 V/cm 
drift field. A 2D Gaussian fit was applied to the S2-S1 profiles to find the anti-correlation axis (black diagonal 
line). Right panel; The combined energy spectrum showing the 511 keV line together with Gaussian fits plus 
linear background are shown to give an energy resolution (CJ/E) of 1.84%. With XYZ position corrections 
applied the energy resolution is 1.62%. 
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Figure 10. Energy resolution for small fiducial volume selection, with XYZ position corrections applied, 
otherwise like Fig. 9. The ultimate energy resolution obtained for the 511 keV line is 1.36%. 


the 511 keV gamma energy line is obtained. It is worth noting that when applying no XYZ position 
corrections, the energy resolution of the 511 keV line is still 1.84%. 

In order to determine the a and j3 values (see Eq. ^71]) , 2D Gaussian fits were applied to the 
S2-S1 profiles of each of the six spectra displayed in Fig. 0 to find their anti-correlation axes. 
Using the BS1 and TS2 intercepts, BS1/ and TS2/, of the anti-correlation axis in each spectrum, a 
and p values were determined for each photo-absorption peak according to a = lF(BSl//£') and 
P = W{TS2j/E), where E is the energy of the photo-absorption events. Averaging over all sources 
yields a = (0.239 ± 0.012) pe/photon and p = (16.1 ±0.6)pe/e^ for the MiX detector at a drift 
field of 200 V/cm. 

The zero field light yield of the detector, when extrapolated from the anti-correlation axis to 
the BSl intercept, is (17.1 ±0.7)pe/keV at 122keV. A correction factor of 0.89 is applied to 
the extrapolated light yield to obtain an estimated true zero field light yield of (15.2 ±0.6) pe/keV 
at 122 keV. This corresponds to 64photons/keV, in agreement with NEST [PIl]. 

The result of high light yield is good energy resolution and discrimination power between 
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Figure 11. Energy resolution in liquid xenon measured by Aprile et al. ||14| |, XENONIOO [|15|| , Lin et 
al. [ P^ and this work. The solid (red) line represents the fit result to the data above 100 keV to a function 
(J/E= (22.2/v^ + 0.42)%. 


nuclear and electron recoils over a large energy range. Energy resolution as a function of energy is 
investigated in this work by making fiducial XY cuts (see Fig. ^ to select events from the central 
part of the detector for each source. The resulting energy spectra for each source were fit with 
background plus gaussian functions. Figure [T^ shows the results from the annular-shaped fiducial 
XY region (bounded by the two red circles) yielding a resolution of a/E = (1.36% ± 0.07)%. 
This region corresponds to a fiducial volume of 1.1 gram of FXe, a small fraction of the entire 
sensitive volume of 114 grams. The restricted annular region is selected for its superior energy 
resolution while still providing enough photo-absorption events to allow a satisfactory signal fit. 
The improvement in energy resolution can be attributed to the axial symmetry of this centralized 
region overcoming the effects of a statistics-limited correction map. 

The energy resolution results for all sources are shown in Fig. Q with an ultimate resolution 
of a/E = (1.03 ±0.11)% at 1.33 MeV, which may be the best achieved in liquid xenon to date. 
Also shown in Fig. m are results obtained in other FXe detectors around the world []I^-|T^. The 
performance of XENONIOO [1T31] is typical of early large scale dual-phase instruments, which did 
not have high photon detection efficiency, while the smaller scale R&D devices represent what can 
be achieved with advancement in TPC design, operation, and sensor performance. The early data 
by Aprile [jHll were collected in a single phase detector using charge amplifiers instead of PMTs 
for the proportional signal. Differences in the resolution results of small-scale R&D devices can be 
attributed to varied quanta collection efficiency from differences in electrode optical transparency, 
TPC aspect ratio, and PMT performance. The excellent energy resolution obtained in the MiX 
detector can be attributed to the very high light collection efficiency, stability in the electric fields, 
PMT gains, and FXe level as well as the ability to reconstruct XY position with high precision, 
allowing for the selection of a fiducial volume that avoids saturation and degeneracies typically 
observed in small FXeTPCs with high symmetry. 
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Table 2. Selection of radioactive sources together with respective energy resolutions are shown. 


Calibration Souree 

Energy 

[keV] 

Resolution 

(%) 

57Co 

122 

2.65 ±0.09 

131mxe 

164 

1.97 ±0.12 

129mxe 

236 

1.70±0.11 

i33Ba 

356 

1.49 ±0.07 

22Na 

511 

1.36 ±0.07 

i3’Cs 

662 

1.36 ±0.06 

^°Co 

1,173 

1.24 ±0.14 

^°Co 

1,333 

1.03±0.11 


4. Conclusions 

We report the design and performanee of the Miehigan Xenon deteetor, a small 3D position sensi¬ 
tive dual-phase liquid xenon time projeetion ehamber with high light yield, long eleetron lifetime, 
exquisite energy resolution over a wide dynamie range and stable operation with minimal human 
interaetion. The ultimate energy resolution of approximately 1% (a/E) at 1.33 MeV may be the 
best resolution aehieved in liquid xenon over a wide energy range to date. If similar energy resolu¬ 
tion is desired in a large deteetor, efforts must be made to improve upon photon detention effieieney 
and inerease sensor eoverage while maintaining eurrent large seale eharge deteetion eapabilities. 

As the Mix deteetor is refined, mueh work on further eharaeterizing LXe as a high fidelity 
radiation imaging medium is on the horizon. A future important goal is to preeisely quantify the 
ER/NR diserimination eapabilities of a EXeTPC under varied eonditions. There are many free 
parameters to optimize in a EXe deteetor, sueh as internal geometry, pressure, drift and extrae- 
tion fields, and EXe level. And fhere is still a quesfion of how fhese paramefers play a role in 
determining fhe ER/NR diseriminafing power in EXe signal produefion. 

Ultimately, fhe modern seienfifie goals of EXeTPCs are fo use fhe feehnology fo diseover 
new physies sueh as WIMP dark mailer and neulrinoless double bela deeay. Physieisls may be 
eloser lhan ever fo aehieving fhese feafs. Additionally, due fo fhe ineredible resolving power of 
fhese defeefors, bolh in energy and in posilion, addilional imporlanl appliealions in EXe-assisled 
tomography and nalional seeurily are wilhin reaeh. 
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Figure 12. Application of the Position-Signal Degeneracy Theorem to the top PMT geometry of the MiX 
detector. Two events occur at points p (in blue) and p' (in red). Loci of points satisfying the equation r = Cr' 
are shown in orange for values of C between 2 and 6. Since the MiX TPMTs (green squares) fall over several 
of these loci, there will be many pairs of indistinguishable points p, p'. 

A. Reconstruction Singularities in Multi-PMT detectors 

It was discovered during analysis of data from a Monte Carlo simulation of the MiX detector that a 
degeneracy exists in the position reconstruction of the detector. Eventually this was formulated as 
a general topological theorem (see Theorem |I]). This degeneracy restricts the ability to reconstruct 
event positions inside a detector in certain regions of space, which in turn may limit the ability 
to make position-dependent signal corrections and fiducial volume selections. Introducing a large 
enough asymmetry into the PMT arrangement can ensure that the degeneracy does not occur. 

A.l Position-Signal Degeneracy Theorem 

For a given event, the number of light quanta that strike each PMT is referred to as the hit pattern. 
Since the hit pattern is determined by the position and energy of the event itself, it is natural to ask 
whether the position and energy of the event are uniquely determined by the hit pattern. In general, 
this is not the case. Even with an arbitrary number of PMTs, it is possible for two events at distinct 
positions to produce indistinguishable normalized hit patterns. If normalization is necessary for the 
position reconstruction model, then theorem |I] applies. A visual representation of how theorem |T] 
limits the position reconstruction capabilities of a detector with MiX-like geometry is illustrated 
schematically in Fig. |T^. 

In the Mix detector geometry, it can be shown that for every S2 event, there is a twin S2 
event that would produce an indistinguishable hit pattern from a different event position. However, 
since the volume in which S2 events may occur is finite, there is a subset of points inside the 
detector whose degenerate partners can never occur. Using numerical methods, one may locate 
these degenerate event positions and place bounds on the spatial region of the detector with good 
position reconstruction. For the MiX detector, the maximum radius with no degeneracy is 29.0 mm 
(as shown in Fig. m- This radius is large enough to avoid position degeneracy in the fiducial 
volume used for fhis analysis (see Fig. 
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Figure 13. Numeric degenerate pairing results from simulated four PMT array. The black curve shows pairs 
of radii r, / in inches from the center of the TPC that are degenerate. Red lines demarcate the maximum 
radius of a MiX S2 event. The green line is r= r' to provide scale. The orange line marks the lowest radius 
for which the degenerate pair radius lies within the MiX S2 region, at 29.0 mm. 


Theorem 1 Let an event occur at some point p G with magnitude M, producing a signal which 
propagates uniformly through space, so that a measuring device with distance r to p will register 
a response of magnitude M/r^. Suppose that the signal magnitude due to this event is measured at 
some points {m, } of distances r,- to p. 

If there exists a point p' whose distances r- to each of the measurement points m,- satisfy rj = Cr, 
for some constant C, then it will be impossible to distinguish an event at point p from an event at 
point p' if the absolute magnitude of the signal is unknown. 


Proof: 

Given the conditions in theorem |^, we may choose magnitude M' = M -C^ such that V/, 


M' M-C^ 






M 


(A.1) 


An event of magnitude M' at the point p' will produce exactly the same signal at each measuring 
point, making it indistinguishable from the event at p. The finite optical position resolution of 
detectors only increases the number of available radii which satisfy this condition. 
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